Measurement of neutrino velocity with the MINOS detectors 
and NuMI neutrino beam 
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The velocity of a ~3 GeV neutrino beam is measured by comparing detection times at the Near 
and Far detectors of the MINOS experiment, separated by 734 km. A total of 473 Far Detector 
neutrino events was used to measure (v — c)/c = 5.1 ± 2.9 x 10 -5 (at 68% C.L.). By correlating the 
measured energies of 258 charged-current neutrino events to their arrival times at the Far Detector, 
a limit is imposed on the neutrino mass of m„ < 50 MeV/c 2 (99% C.L.). 

PACS numbers: 14.60.Lm (Properties of ordinary neutrinos) 
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I. INTRODUCTION 

Investigations of the intrinsic properties of neutrinos 
and their interactions have led to discoveries such as fi- 
nite mass, lepton-flavor number violation, and oscilla- 
tions with large mixing angles These surprises war- 
rant careful measurement of other basic neutrino proper- 
ties such as the relationship between energy and velocity. 

If the mass of the heaviest neutrino is assumed to 
be 3 eV/c 2 , the best direct limit on a neutrino mass 2, 
[|[,then the relativistic velocity v of a 10 GeV neutrino 
should satisfy \v — c\/c < 10 -19 . Cosmological measure- 
ments [|| give a mass limit an order of magnitude smaller, 
implying an even tighter velocity constraint. However, 
theories have been proposed to allow some or all neu- 
trinos to travel along "shortcuts" off the brane through 
large extra dimensions [5|, and thus have apparent ve- 
locities different than the speed of light. Some of these 
theories [f| 0, H[ allow \v — c\/c ~ 10~ 4 at neutrino en- 
ergies of a few GeV. Terrestrial neutrino beams could 
measure an effect of this magnitude. 

Earlier terrestrial measurements [t| [l(| EH constrained 
\v — c\/c < 4 x 1CP 5 by comparing the interaction times 
of muons and muon neutrinos of E v > 30 GeV created in 
a 1 ns beam spill over a distance of ^500 m. This work 
differs in several respects: First, MINOS employs a lower 
energy beam (~3 GeV). Second, we measure the abso- 
lute transit time of an ensemble of neutrinos, to < 100 ns 
accuracy, by comparing neutrino arrival times at the MI- 
NOS Near Detector (ND) and Far Detector (FD), sep- 
arated by a distance of 734 km. Third, we make the 
unique measurement of comparing the energies of neutri- 
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Baseline: 

Distance" ND to FD, L 734 298.6 ±0.7 m [12] 

Nominal time of flight, r 2 449 356 ± 2 ns 

MINOS Timing System: 

GPS Receivers TrueTime model XL-AK 

Antenna fiber delay 1115 ns ND, 5140 ns FD 

Single Event Time Resolution <40 ns 
Random Clock Jitter 100 ns (typical), each site 

Main Injector Parameters: 

Main Injector Cycle Time 2.2 seconds/spill (typical) 

Main Injector Batches/Spill 5 or 6 

Spill Duration 9.7 fis (6 batches) 

Batch Duration 1582 ns 

Gap Between Batches 38 ns 

"Distance between front face of the ND and the center of the FD. 
TABLE I: Relevant MINOS and NuMI Parameters 



nos in charged-current (CC) interactions to the interac- 
tion times in the FD. 

The MINOS detectors QJ, Q HU are steel-scintillator 
tracking calorimeters. Planes of 2.54 cm thick steel sep- 
arate planes made of scintillator strips, 4.1 cm wide and 
1 cm thick. The planes are oriented 3.3° from the normal 
to the beam direction. Strips are aligned orthogonally on 
adjacent planes to allow three-dimensional reconstruc- 
tion of event topology. Multi-anode photo-multiplier 
tubes (PMTs) read out the strips via wavelength-shifting 
optical fibers. 

The time of each PMT hit is recorded by the detec- 
tor's clock to a precision of 18.8 ns (ND) and 1.6 ns 
(FD). Although the implementation of these clocks dif- 
fer, each is synchronized to a Global Positioning Sys- 
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tern (GPS) receiver to provide absolute Universal Coor- 
dinated Time (UTC). The two identical receivers, sit- 
uated underground, connect to the surface by identical 
optical transceivers and optical fibers with lengths given 
in Table [J 

The NuMI beam is created by impinging protons from 
the Fermilab Main Injector onto a graphite target. The 
secondary mesons are focused by two horns and allowed 
to decay in a 675 m long decay pipe. The resulting neu- 
trino beam is 93% v^, 6% V u , and 1% v e + v e at the Near 



Description 



Uncertainty (68% C.L.) 



detector. After oscillating [14j j . the beam at the Far de- 
tector is approximately 60% v^. The energy spectrum is 
peaked at approximately 3 GeV, with a long high-energy 
tail extending to 120 GeV. 

The Main Injector accelerates protons to 120 GeV and 
sends them to NuMI by single-turn extraction. It op- 
erates in one of several modes, allowing either 5 or 6 
batches of protons per spill. A pulsed dipole magnet ex- 
tracts protons from the Main Injector. The extraction 
magnet signal is time-stamped by the ND GPS receiver 
and defines time of the spill, to- 



ll. DATA SELECTION 

CC events in the ND are selected using criteria iden- 
tical to those of Ref. [lij]: events are required to have 
total reconstructed energy less than 30 GeV, have a ver- 
tex contained within aim radius fiducial volume, and 
be in time with the spill (±~7 jits). A probability-based 
particle identification parameter removes neutral-current 
(NC) shower events. This analysis sampled 1.6 x 10 6 ND 
events, roughly 1/3 of the first year's data. 

The pre-selection of events in the FD requires event 
times within ±50 /is of the expected arrival time (as- 
suming a massless neutrino) . Events are accepted if they 
satisfy one of three selections: CC events contained in 
the fiducial volume, neutrino-induced muons from CC in- 
teractions in the rock outside the detector (rock- muons), 
and shower events. 

The contained CC event selection is again similar to 
that of Ref. [T3|; events are required to have a ver- 
tex within the fiducial volume, and to have a well- 
reconstructed track with direction within 53° of that of 
the beam. In this work, both and candidates are 
selected. Events with tracks penetrating the top of the 
detector volume are vetoed as possible cosmic-ray con- 
tamination. The energy of the CC events is determined 
by summing hadronic shower energy and muon energy 
derived from track length or curvature. 

Rock-muon events are selected by considering only 
muons which enter the front face of the detector, to re- 
duce background. The track is required to be contained 
within the detector volume or to exit the lower half of 
the detector, to remove background due to cosmic rays 
reconstructed with the wrong directionality. The track 
was also required to have a direction within 26° of the 
beam direction. 



A Distance between detectors 


2 ns 


B ND Antenna fiber length 


27 ns 


C ND electronics latencies 


32 ns 


D FD Antenna fiber length 


46 ns 


E FD electronics latencies 


3 ns 


F GPS and transceivers 


12 ns 


G Detector readout differences 


9 ns 


Total (Sum in quadrature) 


64 ns 



TABLE II: Sources of uncertainty in v relative time measure- 
ment. 



Shower events are mostly from NC interactions, but 
also include CC events from v ei and v T or in which 
no muons are detected. Events with a cluster of hit 
strips with a total pulse-height greater than approxi- 
mately 300 MeV are accepted as shower events. Shower 
events are required to be contained inside the detector 
volume, thereby reducing cosmic ray events and random 
noise. 

A total of 473 neutrino-induced FD events were se- 
lected, of which 258 were contained Vn or z/„ CC events. 
By relaxing the selection cuts, the cosmic-ray induced 
background is estimated to be < 1 event. 



III. NEUTRINO EVENT TIMING 

The time of a neutrino interaction in the ND is taken 
as time of the earliest scintillator hit, t^D- This time is 
compared to the time of extraction magnet signal, to, and 
corrected for known timing delays: t\ = t^D ~to — d^D- 
Similarly, for FD events, t-i = tpD — to — cLfd- 

The corrections d^D and dFD incorporate known off- 
sets and delays due to readout time, electronic latency, 
and GPS antenna fiber delays. Test-stand measurements 
were used to find the magnitude of each offset. Table ITT1 
summarizes the uncertainties on these corrections. The 
uncertainty on the net correction |d/v£> — <^fd| was deter- 
mined to be ±64 ns at a 68% C.L. The delay of optical 
fibers that run between the surface antennas and the un- 
derground GPS receivers created the largest uncertain- 
ties; these uncertainties were estimated from the disper- 
sion of multiple independent measurements of the fiber 
delays. For example, the delay of a Far Detector fibre 
was independently measured with four instruments: an 
Optical Wavelength Labs BOLT-NL, an Aligent E6000B 
Optical Time Domain Refiectometer (OTDR), an EXFO 
FTB-300 OTDR, and a custom-built test apparatus. The 
results of these four measurements all differed, with an 
RMS of 46 ns, roughly 1% of the delay. 

If the pulse of neutrinos were instantaneous, the devia- 
tion from the expected time-of-flight r could be measured 
as S = (t2 — ti) — t. However, the NuMI beam pulse is 
9.7 (is long, with an intensity time-profile consisting of 6 
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FIG. 1: Neutrino event time distribution measured at the 
MINOS Near Detector. The top plot corresponds to events 
in 5-batch spills P® (ti) while the bottom plot corresponds to 
6-batch spills Pf(ii). 



FIG. 2: Time distribution of FD events relative to prediction 
after fitting the time-of- flight. The top plot shows events 
in 5-batch spills, the bottom 6-batch spills. The normalized 
expectation curves Pi(t) and P£(t) are shown as the solid 
lines. 



batches with short gaps in between. The relative intensi- 
ties of these batches, and the shape of the gaps, are due 
to the proton intensity profile of the NuMI beam. Two 
running modes, '5-batch' and '6-batch' are distinguished 
on a pulse-by-pulse basis. (The '5-batch' mode includes 
two types of spill, where either the first or last batch are 
not delivered to NuMI.) The ND measures this intensity 
profile with neutrino interactions. This measurement cal- 
ibrates the neutrino time relative to the extraction fire 
signal. This profile is shown as Figure [T] and is repre- 
sented as probability density functions (PDFs) Pi{t\) 
and ff(ti). 

The arrival time distribution of neutrinos at the FD 
is similar, but the relative jitter of the two GPS clocks 
further degrades the time resolution. These clocks have 
a maximum error of ±200 ns relative to UTC, with a 
typical error of 100 ns. The uncorrelated jitter of two 
clocks, in addition to detector time resolution, gives a 
total relative (FD/ND) time uncertainty of a — 150 ns. 
We therefore compose PDFs of the expected FD neutrino 
arrival time distribution P% and Pf : 

xP[ l {t')dt' (n = 5,6) (1) 

The resulting PDF describes the predicted time dis- 
tribution at the FD. The time of each event in the FD 
(t\) was compared to this PDF. The deviation 5 from the 
expected time was found by maximizing a log-likelihood 
function (L), summing each event (i) in the 5- and 6- 
batch data: 

L = ^ln P 2 (f 2 -r-5). (2) 



The distribution of measured FD times is shown in 
Figure [H along with the predicted distribution for the 
best fit value of S. The deviation was found to be 
S = —126 ± 32(stat.) ± 64(sys.) ns at a 68% confidence 
limit (C.L.). The systematic uncertainty is due to timing 
offsets shown in Table [TT1 The goodness-of-fit probability 
was determined by Monte Carlo (MC) to be 10%, and 
the likelihood is Gaussian in 5. 



IV. RELATIVISTIC MASS MEASUREMENT 

If the neutrino had a relativistic mass m„ and total 
energy E v , the time of flight would be: 



T mv {E v ) — 




where r is the time of flight of a massless particle. For 
contained z/ M and P p CC events, the MINOS detectors 
measure the neutrino energy, allowing the hypothesis of 
a non-zero mass to be tested by measuring the arrival 
time as a function of E v . 

The measured event times are fitted to a model with 
m v as a free parameter. For symmetry, the fit allowed m v 
to become negative; positive masses indicated retarded 
arrival times, while negative masses were interpreted as 
advanced arrival times, i.e. T- m „ — > t + (r — T m J. The 
detector energy response R{E U , E leco ) was parameterized 
by a PDF derived from MC, where E leco is the energy 
reconstructed in the detector. The true energies El are 
unknown, and so are fitted as a set of 258 parameters 
constrained by R and E rcco . At the peak beam energy, 
R is approximately Gaussian in E v — E rcco with a width 
of -30%. 
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The offset 5 is taken as a parameter constrained by 
the earlier systematic measurements as a Gaussian about 
zero with as = 64 ns. A log-likelihood is constructed 
using the expected arrival time PDF Pzfa), the arrival 
times t\, the fitted true energies E l v . 

L = ^"2 +E t ln P 2 (4 - T mv (El) -6)+ lnR(El,El cco )} 

$ i 

(4) 

The result of the fit was m v = 17±28( stat -) MeV/c 2 at 
a 68% C.L. The likelihood function is non-Gaussian; the 
99% C.L is m u = 17±^(stat.) MeV/c 2 . The best fit gave 
5 = —46 ns and a goodness-of-fit probability of 8%. 

The uncertainty on the energy resolution of the FD 
dominates the systematic uncertainty. The response 
R (El,E* eco ) was found by comparing the reconstructed 
neutrino energy with the input neutrino energy in a 
MC simulation, with events weighted for an oscillated 
neutrino beam, with Am| 3 = 0.0027 eV 2 ,sin 2 20 2 3 
1 .() pjj]. To estimate the systematic uncertainty, the de- 
tector response R was varied by (a) changing the ex- 
pected neutrino energy distribution by varying oscilla- 
tions parameters within the allowed range of Ref. [l4| 
(b) increasing the NC contamination of the CC sample 
by ±50%, (c) changing the shower energy scale by ±11%, 
and (d) changing the muon energy scale by ±2%. For 
each change, R was evaluated and the data re-analyzed. 
We incorporate these systematics by simply taking the 
extremum limits on m v from all of these trials, obtaining 
a final result of m„ = 17i|g(stat.+sys.) MeV/c 2 , 99% 
C.L. The limiting case of m u — 50 MeV/c 2 is shown 
graphically in Figure 02 which shows the data for events 
with energies less than 10 GeV. Neutrinos consistent with 
this mass fall inside the shaded region. 

In practice, this method uses the high-energy events 
to constrain S, and uses the lowest-energy events to con- 
strain the relativistic neutrino mass. If the constraint on 
S is removed, a free fit gives m„ = 14lg(stat.) MeV/c 2 
and 5 = -99 ± 140 (stat.) ns at a 99% C.L., with a 
probability of fit of 10%. 

V. CONCLUSIONS 

By measuring the arrival time of 473 contained CC, 
rock-muon, and NC events as measured by the MINOS 
GPS clocks, the deviation from the expected time at the 
Far Detector was found to be 

5 = -126 ± 32 (stat.) ± 64 (sys.) ns 68% C.L. 

By comparing to the nominal time of flight t, we interpret 
this as a neutrino velocity of v = L/(t + 5) to satisfy 

( v ~ c > = = 5.1±2.9(stat.+sys.) x 10~ 5 68% C.L. 

C T + S V ' 

for neutrinos of ^3 GeV. This measurement is consistent 
with the speed of light to less than 1.8<7. The correspond- 
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FIG. 3: The time and reconstructed energy for contained i/ M 
charged current events. The points show the measured times 
of events and reconstructed energy -B re co • The horizontal error 
bars indicate the ~lcr energy uncertainty. The gray filled 
region indicates the allowed range of times predicted by a 
neutrino with m v = 50 MeV/c 2 . The solid lines indicate the 
allowed range predicted m v — 17 MeV/c 2 . 



ing 99% confidence limit on the speed of the neutrino is 
-2.4 x 10" 5 < (i> - c)/c < 12.6 x 10" 5 . 

This measurement has implicitly assumed that the 
mi and neutrino mass eigenstates that comprise the 
beam are traveling at the same velocity. This assumption 
is borne out in observing that the arrival times at the Far 
detector match the expectation distribution. Indeed, if 
the two eigenstates were to travel at velocities differing 
by as little as Av/v >4x 10 -7 , the short ^1 ns bunches 
would separate in transit and thus decohere, changing or 
destroying oscillation effects at the Far detector. 

Besides the novelty of the technique, this measurement 
is unique in that it probes the 1-30 GeV region of neutrino 
energy not measured by previous experiments. The mea- 
surements described in Refs. [1, reached a sensitivity 
slightly better than this work, but only for neutrinos of 
~ 25 GeV and higher. The most sensitive test of neutrino 
velocity was achieved by comparing 1 the arrival times 
of neutrinos [ill 13 and photons from SN1987a :i which 
achieved a sensitivity of \v — c\/c < 2 x 10~ 9 [H, [20| . 
four orders of magnitude better than the terrestrial mea- 
surements, but only for neutrinos of energy ~ 10 MeV. 
In principle, neutrino velocity could be a strong function 
of energy. Our measurement constrains this previously 
untested regime. 

By using the arrival time and reconstructed energies 
of 258 contained CC and events, a limit was found 



The SN1987A measurement is predicated on the theoretical as- 
sumption that neutrinos and photons are emitted within three 
hours of each other. 
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on the relativistic mass of the neutrino: 

m v < 50 MeV/c 2 (stat.+sys.) 99% C.L. 

The method of relating time-of-fiight to neutrino en- 
ergy is a new technique made available by the MINOS 
calorimeters. With the entire MINOS data sample we 
anticipate a factor of 10 increase in statistics. This will 
improve the sensitivity of this measurement to as low as 
10 MeV/c 2 . 

Direct mass measurement yields much tighter con- 
straint of neutrino mass 0, Q , but the method of mea- 
suring mass by time-of-fiight can additionally be viewed 
as a test of the relativistic energy-velocity relation. If 
the relation holds, the mass measured by time-of-fiight 
should be consistent with direct measurements. Alter- 
nate theories, such as those suggested by Refs. [f| 0, H| 
could be tested by these data. 
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